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Abstract 

The observations of the surfaces of the mid-sized Saturnian satellites made 
by Cassini-Huygens mission have shown a variety of features that allows 
study of the processes that took place and are taking place on those worlds. 
Research of the Saturnian satellite surfaces has clear implications not only 
for Saturn’s history and Saturn’s surroundings, but also for the Solar System. 
Crater counting from high definition images is very important and could serve 
for the determination of the age of the surfaces. In a recent paper, we have 
calculated the production of craters on the mid-sized Saturnian satellites by 
Centaur objects considering the current configuration of the Solar System. 
Also, we have compared our results with crater counts from Cassini images by 
other authors and we have noted that the number of observed small craters 
is less than our calculated theoretical number. In this paper we estimate 
the age of the surface for each observed terrain on each mid-sized satellite 
of Saturn. All the surfaces analyzed appear to be old with the exception 
of Enceladus. However, we have noticed that since there are less observed 
small craters than calculated (except on lapetus), this results in younger ages 
than expected. This could be the result of efficient endogenous or exogenous 
process(es) for erasing small craters and/or crater saturation at those sizes. 
The size limit from which the observed number of smaller craters is less than 
the calculated is different for each satellite, possibly indicating processes 
that are unique to each, but other potential common explanations for this 
paucity of small craters would be crater saturation and/or deposition of E- 
ring particles. These processes are also suggested by the Endings that the 
smaller craters are being preferentially removed, and the erasure process is 
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gradual. 

On Enceladus, only mid and high latitude plains have remnants of old 
terrains; the other regions could be young. In particular, the regions near 
the South Polar Terrain could be as young as 50 Myrs old. On the contrary 
for lapetus, all the surface is old and it notably registers a primordial source 
of craters. As the crater size is decreased, it would be perceived to approach 
saturation until D <2 km-craters, where saturation is complete. 
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1. Introduction 

The surfaces of mid-sized Saturnian satellites are a kind of laboratory 
where it is possible to observe and investigate the physical and dynamic 
processes that have been taking place around Saturn and also throughout 
the Solar System. The study of impact craters on the satellite surfaces allows 
us to better understand what could be happening in the Saturn environment 
to produce these craters. The Saturnian satellite system was observed and 
studied in the past by Pioneer 11 and Voyager 1 and 2 spacecrafts, greatly 
increasing the knowledge of the Saturn system. Voyager images revealed 
diverse satellite surfaces (Smith et ah, 1981, 1982). Crater counts on those 
images indicated that Saturnian satellites have been variably cratered, which 
suggest different geologic histories (Plescia and Boyce 1982, 1983, 1985). At 
present, the Cassini-Huygens mission is visiting the Saturn system, and the 
detailed observations it renders provide us with new paradigms and physical 
processes to understand and interpret. 

The mid-sized icy satellites of Saturn are Mimas, Enceladus, Tethys, 
Dione, Rhea, and lapetus. They are regular satellites, mainly composed 
of water ice and in synchronons rotation. The Cassini-Hnygens mission has 
observed all of them in detail allowing scientists to obtain accnrate infor¬ 
mation on the shapes, mean radii and densities (Thomas, 2010) and gravity 
fields (Jacobson et ah, 2006). Fnrthermore, some of these satellites show 
traces of physical activity and renovation, possibly dne to recent geologi¬ 
cal processes. The geologic activity could include endogenous activity such 
as viscous relaxation, volcanism, and/or tectonic or even atmospheric pro¬ 
cesses. Cratering itself is a potential process in which the formation of large 
craters could remove small craters by ejecta emplacement and seismic shak¬ 
ing. Also, exogenous processes, such as in fall from the E-ring or debris rings. 
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could erase surface features. 

Those physical processes would renew the satellite surfaces, erasing old 
or young surface features such as craters. The most striking case is Ence- 
ladus whose surface differs markedly from region to region and has present 
active geysers emanating from four parallel fractures in the south polar re¬ 
gion, called “tiger stripes” (Porco et ah, 2006). This region is a very active 
one where the surface is modihed, erasing the craters (totally or partially). 
On Dione, in turn, Cassini observations reveal two different surfaces: the 
cratered plains (cp), heavily cratered, and the smoother plains (sp), with 
a lower cratering suggesting a younger surface. On Rhea, measurements by 
Cassini spacecraft detected a tenuous atmosphere of oxygen and carbon diox¬ 
ide (Teolis et ah, 2010). lapetus is the opposite case of Enceladus as it has 
heavily cratered plains with large degraded basins, which indicates that its 
surface is ancient. What is more, the count of small craters showed a size 
distribution indicative of crater saturation (Denk et ah, 2010). 

Kirchoff and Schenk (2009) and (2010), hereafter KS09 and KSIO, an¬ 
alyzed high-resolution Cassini images and obtained the number and size- 
frequency distribution (SFD) of craters for the mid-sized icy satellites. With 
their observed number of craters and the previous cratering rate estimations 
by Zahnle et al. (2003), they calculated the surface ages of each satellite for 
some crater diameters. 

There are a number of factors that must be taken into account in the 
analysis of the age of a satellite area. First, all possible impactor popula¬ 
tions should be considered. The analysis of the images obtained by Voyager 
(Smith et ah, 1981, 1982) implied that the satellites were struck by two dif¬ 
ferent impactor populations: Population I, which produced a greater number 
of large craters (bigger than 20 km), and Population II, which produced a 
greater number of smaller craters (smaller than 20 km) (Smith et ah, 1982). 
The origin or even the existence of both populations is disputed. Smith et al. 
(1981, 1982) suggested that Population I was the tail-off of a postaccretional 
heavy bombardment, while Population II has the form expected for collisional 
debris from the satellites or other orbiting debris. Horedt and Neukum (1984) 
concluded that cratering on Saturnian satellites is produced by heliocentric 
objects as well as by planetocentric impactors. On the other hand, Hartman 
(1984) noted that crater densities on heavily cratered surfaces throughout 
the Solar System are all similar due to a “saturation equilibrium”. He ar¬ 
gued that only one single population of heliocentric planetesimals and their 
fragments had been recorded. Dobrovolskis and Lissauer (2004) studied the 
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fate of ejecta from the irregularly shaped satellite, Hyperion, suggesting that 
it does contribute to Population II craters on the inner satellites of Saturn. 
However, those particles would produce craters with a different morphology 
than those produced by a heliocentric source. 

Moreover, consideration of the primordial sitnation is important for the 
study of the origin of craters. It is believed that the mass of the primor¬ 
dial trans-Neptnnian zone, the source of Centaurs, was ~ 100 times higher 
than the present one (Morbidelli et ah, 2008). This mass might have been 
depleted by a strong dynamical excitation of the trans-Neptnnian region. 
There were several models that described the mass depletion; for example, 
in the “Nice Model” and its snbseqnent versions, the interaction between the 
migrating planets and planetesimals destabilized the planetesimal disk and 
scattered the planetesimals all over the Solar System before the time of the 
LHB (Tsiganis et ah, 2005; Levison et ah, 2008). Accordingly, there was 
mnch primordial mass that struck the planets and their satellites very early 
in the Solar System history. It wonld be expected that this event marked the 
surfaces of the satellites, mainly producing a great nnmber of larger larger 
craters than at present. Those craters are probably the larger ones that can 
be observed on the satellites. Moreover, there are papers (Nervorny et ah, 
2003, 2007) that argue in favor of more primordial irregnlar satellites and a 
more active primordial collisional activity aronnd the major planets. There 
would even be populations of large irregnlar satellites and debris of catas¬ 
trophically disrnpted satellites that may have played important roles in the 
history of the Saturn system (Dones et ah, 2007). 

A primordial crater contribntion to the Satnrn system could also be con¬ 
nected with the crater satnration observed on some of the Satnrnian satellites. 
The formation of craters, especially the larger ones, is a process that is effec¬ 
tive in erasing small craters, not only through crater formation, bnt also by 
the ejecta blanket and seismic shaking. When a snrface is so heavily bom¬ 
barded that the formation of craters is eqnaled by the obliteration of craters, 
it reaches a crater satnration equilibrinm. In this case, the density of craters 
on a surface does not change proportionally. This process is then critical to 
the determination of the sonrce of impactors, geological processes, ages, etc 
(Hartmann and Gaskell, 1997). Crater saturation eqnilibrinm hnally causes 
a distinctive cumnlative size freqnency distribntion (SFD) of craters. Gault 
(1970) obtained a —2 power-law, from a model based on first principles and 
laboratory experiments. Hartmann (1984) htted a —1.83 power-law for the 
cnmulative SFD to crater connt data of the snrface of various Solar System 
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objects. Richardson (2009) developed a detailed cratered terrain evolution 
model and observed the way in which crater densities attain equilibrium con¬ 
ditions. He found that if the impactor population has a cumulative power-law 
slope of < —2, crater densities reach a cumulative power-law of about —2; and 
if the impactor population has a cumulative power-law slope of > —2, crater 
density equilibrium values follow the shape of the production population. 
Squyres et ah (1997) argued that cratering on a surface is a random process 
and as crater obliteration becomes more important, the spatial distribution 
of craters tends to be uniform. They used both the spatial distribution and 
size-frequency distribution of craters to study, with statistical techniques, 
how a cratered surface approaches saturation. They found that at least 25% 
of the craters on Rhea and Callisto were destroyed by subsequent oblitera¬ 
tion. Therefore, crater saturation is especially important on highly cratered 
surfaces that are in general old. 

In the supposedly young surfaces, as in the case of Enceladus for exam¬ 
ple, geological processes are erasing craters. However, such processes compete 
with possible exogenous processes, like particle deposition on the surface that 
can £11 in the craters, finally removing them. Mid-sized Saturnian satellites 
with the exception of lapetus are embedded in the Saturn E-ring; therefore 
their particles cross the orbits of the satellites with some probability of im¬ 
pact. Enceladus’s plumes are thought to be the source and maintenance of 
E-ring material (Hamilton and Burns, 1994, Kempf et al., 2010). A frac¬ 
tion of the plume particles escape to populate the E-ring but other fraction 
returns to Enceladus hitting its surface (Kempf et al., 2010). Ingersoll and 
Ewald, (2011) estimated (12 ± 5.5)10® kg for the mass of particles in the 
E-ring and a particle lifetime in the E-ring of about 8 years. This mate¬ 
rial could be an exogenous source causing erasure of craters in the mid-sized 
Saturnian satellites. 

In a recent paper, Di Sisto and Zanardi (2013), hereafter DZ13, calculated 
the production of craters on the mid-sized Saturnian satellites produced by 
current Centaur objects coming from the Scattered Disk (SD) and plutinos, 
in the trans-Neptunian region, and compared these calculations with the 
Cassini observations made by KS09 and KSIO. Also obtained was the current 
cratering rate on each satellite. In that paper, the authors concluded that 
since the number of observed small craters is lower than their calculated 
theoretical number, and because there are physical processes on at least 
some satellite surfaces that erode them, such difference is likely to be caused 
by those processes. It seems that the satellite surfaces were “reset” by those 
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physical processes and were, therefore, younger and dehned by the current 
geological processes. However, as mentioned, a number of factors has to be 
considered to test what processes are really taking place on the satellites and 
how these affect the age of their surfaces. Therefore, in this paper we analyze 
in detail our theoretical predictions about the current cratering production 
to obtain satellite surface ages. These age determinations are not only a hrst 
approximation but also a way to help determine the time scale of the dynamic 
and geological processes that are taking place on the mid-sized satellites of 
Saturn. 

Since each satellite is a unique world, it is also interesting to test pecu¬ 
liarities and different processes on each satellite in our study. 

2. Age calculation 

For the calculation of the number of impact craters and ages on the 
mid-sized Saturnian satellites all potential populations of impactors should 
be considered. Those impactors could be main-belt asteroids, Jupiter and 
Neptune Trojans, Jupiter Family Comets (JFCs) and Centaurs from the 
SD, irregular satellites, planetocentric bodies and Nearly Isotropic Comets 
(NICs) (Dones et al. 2009). There is evidence, such as recent collisions and 
passages of comets by Jupiter and the observations of Saturn-crossing comets 
and Centaurs (Dones et ah, 2009) as well as also previous work by Zahnle et 
al. 1998 and 2003, that Centaurs and JFCs dominate impact cratering in the 
Outer Solar System. The only truly unknown potential impactor population 
is the planetocentric one. This hypothetical source has been discussed by 
previous papers (see Sect. 1) but it has not been observed so far; its possible 
contribution or existence will be analyzed for each satellite with our results. 

Di Sisto and Zanardi (2013) calculated the production of craters on the 
mid-sized Saturnian satellites by Centaur objects from SD and plutinos. 
They found that the contribution of plutinos is negligible with respect to 
Scattered Disk Objects (SDOs). They used a method previously developed 
by Di Sisto and Brunini (2011), hereafter DBll. Both papers are based on 
the numerical simulation carried out by Di Sisto and Brunini (2007), who 
studied the evolution of SDOs for 4.5 Gyrs when they enter the Centaur 
zone, considering the current configuration of the Solar System. Di Sisto 
and Brunini (2007) studied only Centaur objects when they enter the giant 
planet zone from the SD; they stopped their simulation when Centaurs enter 
the JFCs zone. Thus, the contribution of JFCs to the satellite cratering is 
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not considered here. However, only a ~ 22% of SDOs enter the JFC zone 
(Di Sisto and Brunini, 2007). These objects were studied by Di Sisto et ah 
(2009) through a numerical simulation that included a complete dynamical- 
physical model. Only a fraction of the initial comets return to the Centaur 
zone, depending on the size of the comet (a JFC reduces its radius due to 
sublimation and splitting). We can then assume that this contribution should 
be approximately one order of magnitude smaller than that of the Centaurs 
when they enter the Saturn zone from the SD. 

DZ13 and DBll used the output files of the encounters of SDOs with 
Saturn when they enter the Centaur zone (from Di Sisto and Brunini, 2007) 
to calculate the cumulative number of craters (iVc(> D)) produced by current 
Centaurs on the Saturnian satellites. It is generally accepted that the mass 
of the primordial trans-Neptunian zone, the source of centaurs, was ~ 100 
times higher than the present one (Morbidelli et ah 2008 ) and decreased to 
its present value in at most 1 Gyr (“Nice Model”). This early phase of the 
Solar System was not modeled in the numerical simulation by Di Sisto and 
Brunini (2007), which instead started after the initial mass depletion, when 
the Solar System began to stabilize, and essentially becomes like the current 
conhguration of the Solar System. As mentioned, this simulation was run for 
4.5 Gyrs, since the exact duration of this early phase is not known; but it 
should be representative of the post initial mass depletion time, namely < 4 
Gyr. Then, Nc{> D) corresponds to the number of craters produced since 
the change in conhguration of the Solar System predicted by the Nice Model 
(~ 4 Ga). 

In DZ13, we calculated the W(> D) for each satellite and for two cases 
of the SFD of SDOs. That is, we considered a SFD of SDOs as a power law 
with a differential index si = 4.7 for diameters d > 60 km (from Elliot et ah, 
2005). The SFD breaks for d < 60 km, but given the uncertainty in the SFD 
for small objects, we considered two values of the differential index: S 2 = 2.5 
and S 2 = 3.5 (the Dohnanyi exponent). In DZ13, we also compared our 
theoretical number of craters with the number of observed craters obtained 
by KS09 in Figs. 2-7 of their paper. We noticed that our calculated number 
of craters for S 2 = 3.5 is in general closer to the observed number by KS09. 
However, for small craters, the observed number of craters is lower than the 
calculated one. Since in general there would be physical processes on mid¬ 
sized Saturnian satellite surfaces that erode them, these processes could erase 
small craters. The extent of this erasure is specihc to each satellite because it 
depends on the process that is taking place on each one. It is possible, then. 
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to account for this process by calculating the age of each satellite surface. 
The calculation of the age of each terrain is as follows: hrst, let us assume 
(on the basis of previous comments)that the Centaurs from the SD are the 
main source of craters in the current conhguration of the Solar System and 
that, given the similarity with observations mentioned before, the calculated 
number of craters is the one obtained from the differential index: S 2 = 3.5 of 
the SFD of small SDOs. This assumption can be considered a hrst step in the 
determination of ages. Since Centaurs come mainly from the SD (Di Sisto and 
Brunini, 2007), they should be the most important current heliocentric source 
of craters. It should be highlighted that a possible planetocentric source of 
craters has been excluded since there are no estimations or calculations from 
this hypothetical source. 

From DZ13, we know the total number of craters on each satellite greater 
than a given diameter D, produced by Centaurs from the SD in the current 
conhguration of the Solar System {Nc{> D)). We took the observed number 
of craters W(> D) from KS09. Then, if W(> D) > Nc{> D), we can 
state that the terrain has preserved craters that were produced at earlier 
times of the Solar System, or that there have been other important sources of 
craters (such as a planetocentric source), and/or there is no current geological 
activity that could erode the craters. But if No{> D) < Nc{> D), there are 
Nc{> D) — No{> D) craters that have been erased and then there must be 
a geological process taking place on the satellite surface. The comparison 
between the calculated and the observed numbers can be seen in Figs. 2- 
7 in DZ13, where it is noticeable that there is a range of crater diameters 
(different for each satellite) for which No{> D) < Nc{> D). Therefore, those 
terrains could have a geological activity that has eroded the craters, and the 
surface would be young. 

To calculate the age of those surfaces, the temporal dependence of crater¬ 
ing was used. According to the method developed by DBll, the cumulative 
number of craters produced by Centaurs is proportional to the number of 
encounters of Centaurs with Saturn. Then, the temporal dependence of cra¬ 
tering is equal to that of encounters with Saturn (see Eqs. 2-7 in DBll). 
In Fig. [T]of this paper, the fraction of encounters with Saturn (number of 
encounters for a given time with respect to the total number of encounters) 
was plotted as a function of time and found that the whole plot could be 
htted by a log-function given by: 


F{t) = a logt + b, 


( 1 ) 



t [Myrs] 

Figure 1: Fraction of encounters of SDOs with Saturn versus time and the fit to the data. 

where a = 0.198406 ± 0.0002257 and b = -3.41872 ± 0.004477. 

Then, the cumulative number of craters on a given satellite depending on 
time can be obtained from: 

N,{> D,t) = F{t)N,{> D), (2) 

Differentiating this equation with respect to time, we obtain the cratering 
rate over time from: 


C(,> = D), (3) 

Zanhle et al. (1998) calculated cratering rates on the Galilean satellites 
considering the contribution of several heliocentric sources. They also cal¬ 
culated the age of the satellite surfaces. To such end, they considered that 
the cratering rate on each satellite decreases with time as t~^ from a numer¬ 
ical simulation developed by Holman and Wisdom (1993) who found that a 
collisionless Kuiper Belt dissipates as t~^. This is the same behavior of our 
calculated cratering rate of Eq. [3l since the original simulation by Di Sisto 
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and Brunini (2007) deals with the dynamical evolntion of TNOs. Then, we 
follow a development analogous to Zanhle et al. (1998) to hnd the age of a 
satellite surface but now considering Eq. ([3]) for the cratering rate. Thus, the 
estimated cratering timescale or the “age” of the surface r would be given 
by: 

No(>D) 

t{> D) = tQ{l — e «JVc(>D) (^ 4 ^ 

where to = 4.5 Gyrs is the age of the Solar System. 

r represents the model age of the surface for a given diameter, because it 
is a measure of the cratering timescale, taking into account the production 
of craters and the erosion of the terrain. 

In the following sections, we show our results for the mid-sized Saturnian 
satellites. 

3. Results 

As mentioned in the previous section, considering that Centaurs are the 
current main source of craters, if No{> D) < Nc{> D), there are Nc{> 
D) — No{> D) craters that could have been erased. This difference can be 
plotted against D to get an understanding of the magnitude of the process on 
each satellite. In Fig. [2l we plot Nc{> D) — Ao(> D) versus D for Mimas, 
Tethys, Dione and Rhea considering No{> D) from KS09. For all those 
satellites, the smaller D is, the greater the difference between observation 
and our model. Notice also the different slopes of the curves. Tethys and 
Mimas have more flattened curves meanwhile the curves for Dione and Rhea 
are steeper. Therefore, for Tethys and Mimas, there seems to be a process 
that is acting for a wider range of diameters than for Rhea and Dione. 

We calculate the age depending on D from Eq. [5] and plot the results 
for Mimas, Tethys, Dione and Rhea in Fig. [21 These age-curves correspond 
to the terrains analyzed by KS09. In all the curves, there seems to be a 
noticeable correlation between older age and diameter; that is, smaller craters 
are younger than greater ones. As can be seen, only small craters would be 
young or, at least, younger than 4.5 Gyrs old. This means that there could 
be a geological or other physical process that gradually erases craters and 
modihes the satellite surface; naturally, smaller craters are cleared before 
greater ones. Furthermore, the slopes of the age curves are different for each 
satellite. This fact could imply either that each satellite presents different 
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Figure 2: Difference between the cumulative number of theoretical craters and observed 
craters for Tethys, Dione, Mimas and Rhea according to D. 

geological processes or exogenous erasing agents (like deposition of E-ring 
particles), or at least, that each process has a different temporal scale. It 
could also be possible that on a given satellite surface the equilibrium crater 
saturation had been reached. 

There also exists a correlation between the size range or erased craters 
and distance to Saturn. In the ranges of the diameters plotted, Mimas, the 
innermost mid-sized satellite located in the inner edge of the E-ring, has a 
wider range of diameters of erased craters than Tethys, which is farther away 
from Saturn. Then follows Dione in distance to Saturn, and hnally Rhea, at 
the end of the E-ring. This correlation is very interesting, and potentially 
implies that one can associate a modihcation surface agent that depends on 
the distance to the planet. 

Using Eq. 0 there is a certain diameter D from which No{> D) gets 

No{>D) 

greater than D) and then e “^cOd) approaches 0. Therefore, the cal¬ 

culated ages trend asymptotically to the age of the Solar System and for 
those diameters, the surface would be old. As mentioned in the previous 
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Figure 3: Age of the surface of Tethys, Dione, Mimas and Rhea according to D. 

section, when No{> D) > Nc{> D), there are more craters on the surface 
of the satellite than our calculated number. Since Nc{> D) represents the 
crater contribution from Centaurs in the current configuration of the Solar 
System, there are craters that were produced earlier in the Solar System 
history and/or there is another source of craters, such as planetocentric ob¬ 
jects. The limit diameters beyond which No{> D) > Nc{> D) (or before 
which No{> D) < Nc{> D)) are presented in table [U it can be seen that 
those limit diameters are different for each satellite and might be associated 
with different kinds of processes. Therefore, we cannot evaluate other crater 
sources for all satellites in general, but we can analyze all our results for each 
satellite in separate subsections. 

3.1. lapetus 

lapetus is Saturn’s third-largest moon after Titan and Rhea. Its orbit 
around Saturn has a semimajor axis of 3.56 x 10® km (~ 59 Saturn radii), low 
eccentricity and a slightly high inclination of 8.3°. Cassini images revealed 
a near equatorial ridge system that extends for more than 110° in longitude 
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Table 1: Limit diameter beyond which No{> D) > Nc{> D)) 


Satellite 

D [km 

Mimas 

15 

Enceladus 

none 

Tethys 

25 

Dione-cp 

5 

Dione-sp 

17 

Rhea 

5 

lapetus 

2 


and that rises more than 20 km above the surrounding plains (Porco et 
ah, 2005). The ridge morphology is consistent with an endogenous origin 
(Ip, 2006). Giese et al. (2008) found substantial topography on lapetus’ 
leading side from Cassini images with heights in the range of —10 km to 
+ 13 km. A distinctive feature of lapetus is its global albedo dichotomy. The 
leading side is dark with albedo values of ~ 0.04 on a roughly elliptical area 
(named Cassini Regio after the astronomer Jean Cassini who first noticed 
it); the trailing side, in turn, is relatively bright with albedos of ~ 0.6. 
Two classes of theories have been suggested to account for the origin of this 
dichotomy: the endogenous ones (Smith et ah, 1981, 1982) and the exogenous 
theories which suggest that debris or dust material from the external moons 
of Saturn impact the dark hemisphere (Cruikshank et ah, 1983; Buratti and 
Mosher, 1995). A model including both endogenous and exogenous causes 
was developed by Spence and Denk (2010) which demonstrates that the 
dichotomy can be explained by runaway global thermal migration of water 
ice, triggered by the deposition of dark material on the leading hemisphere. 
From a dynamical point of view, the routes of possible dust impactors on 
lapetus were studied by Leiva and Briozzo (2013). They analyzed low-energy 
incoming dust particles and obtained their distribution on the surface of 
lapetus. There is also a global color dichotomy found on lapetus (Denk 
et ah, 2010), the leading side being redder than the trailing side. Spencer 
and Denk (2010) provided evidence for an exogenous origin for the redder 
leading-side parts and suggested that it initiated the thermal formation of 
the global albedo dichotomy. 

lapetus is heavily cratered, with degraded basins with diameters nearly 
the size of lapetus’s radius, which indicate that both the bright and the 
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dark areas are ancient. Plescia and Boyce (1983) analyzed Voyager 2 images, 
determining the number of craters at large diameters (because of low image 
resolution) that indicated that the bright terrain is an ancient surface that 
dates to the LHB (Plescia and Boyce, 1985). Denk et al. (2010) plotted the 
cumulative crater size-frequency distribution of lapetus No{> D), combining 
hve individual measurements (from Cassini images) that we reproduce in 
Fig. m (extrapolating their counts to the complete lapetus surface) together 
with Plescia and Boyce (1983) counts. In this figure, we also plotted the 
cumulative number of craters on lapetus greater than a given diameter D, 
produced by current Centaurs from the SD (iVc(> D)) calculated from our 
model. As mentioned, we assumed a SFD of SDOs that breaks for d < 60 km 
considering two values of the differential index: S 2 = 2.5 and S 2 = 3.5. The 
two resulting curves are plotted in Fig. 01 As noticed by DZ13 for the other 
satellites, our theoretical number of craters for S 2 = 3.5 (green curve) is closer 
to the observed number. Moreover, both Nc{> D) and No{> D) are similar 
for H < 2 km and for H > 2 km No{> D) > Nc{> D). We already noticed in 
DZ13, comparing our results with those of KSIO, that No{> 5 km) > Nc{> 5 
km), which is consistent with an old surface. Denk et al. (2010) found that 
for D < 5 — 10 km, the crater cumulative SFD follows a —2 power law that 
corresponds to crater equilibrium saturation. In fact, by fitting power-laws 
of the form No{> D) = cD^ to the observed crater distributions plotted in 
Fig. m we found for each terrain and ranges of diameters, different values of 
b. These are shown in Table [21 

In the anti-Saturn Hemisphere region, lapetus-dark and even large basins, 
there seems to be a tendency of a steeper slope for greater sizes. There¬ 
fore, considering that for D < 2 km, the crater distribution is saturated, 
a decrease of b for larger craters might suggest that larger craters are not 
likely saturated. If this interpretation were correct, then the crater distri¬ 
bution would reflect the impactor distribution only for large craters. Our 
theoretical crater distribution Nc{> D) with S 2 = 3.5 (an impactor popula¬ 
tion with a cumulative SFD with an index of 2.5) has a power-law index of 
b = —2.976 for D > 0.08 km. Hence, only for dark terrain and D > 10 km 
of the anti-Saturn hemisphere, the observed crater distribution approaches 
our theoretical crater distribution. On the bright hemisphere, however, the 
observed crater power-law index is —1.7 from Cassini counts and —2.02 from 
Voyager counts, a much flatter slope than the theoretical law. 

On the other hand, from Fig. 01 we note that for D > 5 km, No{> 
D) > Nc{> D) with a difference ~ 10 to ~ 30 times greater. Besides, ac- 
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Table 2: Power-law indexes of the observed cumulative crater distribution obtained by 
fitting to the observed craters by Denk et al. (2010), Plescia and Boyce (1983) and KSIO 


Region 

D-range [km 

b 

Anti-Saturn hemisphere 

5 < D < 35 

-2.2 


8 < D < 35 

-2.36 


10 < D < 35 

-2.47 

Center of basin Engelier 

2 < D < 14 

-1.66 

Large craters on trailing side 

17 < D < 202 

-1.7 

Basins (global) 

70 < D < 350 

-0.92 


350 < D < 684 

-3.48 

lapetus-dark (KSIO) 

4 < D < 80 

-2.66 

lapetus-bright (KSIO) 

4 < D < 65 

-1.7 

lapetus-bright (PB83) 

20 < D < 142 

-2.02 


cording to our calculations, the greatest crater produced by current centaurs 
has D ~ 100 km, but there is observed a signihcant number of craters and 
basins greater than this diameter. Therefore, since our model represents the 
present contribution to cratering, and a current planetocentric population of 
big objects near lapetus is not expected, the main crater source of lapetus 
must be primordial, which is consistent with previous studies (Plescia and 
Boyce, 1983, 1985). As mentioned, it is believed that the mass of the pri¬ 
mordial trans-Neptunian zone was ~ 100 times higher than the present one 
(Morbidelli et ah, 2008), in agreement with the differences found between 
observation and present modeled contribution. 

Therefore, from the comparisons of our model and observations, we can 
conclude that big craters on lapetus are old and connected with the primor¬ 
dial epoch of the Solar System. As the crater size decreases, it seems to 
approach saturation until H < 2 km-craters, where saturation is complete 
(Denk et ah, 2010). This crater saturation does not allow us to evaluate 
the possibility of another current source of craters such as planetocentric 
small objects. The different observed power-laws between bright and dark 
terrains, which may be an artifact of analyzing small areas (Kirchoff and 
Schenk, 2010), cannot be explained by our model. If real, this difference 
could be connected to primordial times or also to probable apex-antapex 
asymmetry (Zhanle et ah, 2003). 
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Figure 4: Cumulative number of craters of lapetus with respect to D. Points represent the 
number of craters obtained by Denk et al. (2010) and Plescia and Boyce (1983) (the gray 
ones), the red curve represents the cumulative number of craters obtained by our model 
for S 2 = 2.5, and the green one for S 2 = 3.5. 


Furthermore, an interesting result obtained from our model is the recent 
production of small craters. Denk et al. (2010) identified brighter-than- 
average craters on highest-resolution images on the dark side, with diameters 
of up to ~ 200 m, the brightest one being the crater Escremiz, which is ~ 4 
times brighter than its surroundings. They estimated crater ages assuming 
that brightness and age are correlated and that the darkening process acts 
uniformly for all fresh craters. They also stated that a new crater is ~ 10 
times brighter than the surrounding dark terrain and that it declines to half 
brightness in ~ 10000 years. Therefore, crater Escremiz is estimated to be 
~ 10000 years old. Moreover, they estimated that on any of lapetus’ dark 
hemisphere, the largest crater with an age similar to Escremiz should have 
D ~ 200 m, and that slightly more than 100 craters of 60 m or larger and 
with an age similar to or younger than Escremiz should exist. From our 
model, it is possible to calculate the time the production of craters with a 
diameter greater than a given D. From Eqs. |2] and [Hand considering that 
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Cassini Regio covers ~ 40% of lapetus surface, we found that the last crater 
greater than D was produced At years ago given by: 

At(> D) = to(l - (5) 

Then, the last crater with D > 60 m in Cassini Regio was produced 
~ 30 years ago, and the last 100 craters with H > 60 m in Cassini Regio 
were produced in the last ~ 3000 years. The last crater with D > 200 m 
in Cassini Regio was produced ~ 800 years ago. Those ages are somewhat 
smaller than the Denk et al.’s (2010) predictions. This might suggest that 
our calculations of the contribution of Centaurs to the recent cratering are 
more numerous than expected. However, these differences could be related 
to the probable apex-antapex asymmetry (Zhanle et ah, 2003), which is not 
considered in this model, or to the intrinsic errors of the current number 
of SDOs at small sizes or even to errors in Denk et al.’s (2010) model of 
age-brightness correlation. 

3.2. Mimas 

Mimas is the innermost mid-sized Saturnian satellite and its orbit is lo¬ 
cated within the inner edge of the E-ring. Mimas has a diameter of 394.4 
km, and has a density of 1.149gr/cm^ slightly higher than water. Mimas’ 
surface is fully covered with impact craters and is one of the most heavily 
cratered bodies in the Solar System. This satellite has a large impact crater 
of 139 km in diameter (one third the diameter of Mimas) called Herschel. 
Schmedemann and Neukum (2011) also suggest evidence for two highly de¬ 
graded large craters, one of which is 153 km in diameter. After calculating 
the production of craters by Centaur objects on Mimas, DZ13 found that the 
largest crater produced by Centaurs has a diameter of 113 km, whose size 
is similar to Herschel crater. Lastly, KSIO used high resolution images from 
Cassini to record surface impact craters of the surface and they found a high 
density of craters with diameters ~ 10 < D < 30 km. Plescia and Boyce 
(1982) obtained crater counts in different areas of Mimas from Voyager im¬ 
ages. They distinguished an equatorial area near Hershel (of low resolution), 
and south polar and adjacent areas which generally lack craters > 20 km. 

Buratti et al. (2011) have failed to hnd plume activity on the satellite 
and concluded that there could only be low-level geologic activity. They also 
limited this plume activity to an upper level of one order of magnitude less 
than the production for Enceladus. Measurements of Mimas’s forced libra- 
tions by Tajeddine et al. (2014) conhrms the libration amplitudes calculated 
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Figure 5: Age of the south polar (PB-SP) and intermediate latitudes regions (PB-T) 
from Voyager and Cassini images (KS09) of Mimas with respect to D. 


from the orbital dynamics, with the exception of one amplitnde that depends 
on Mima’s internal structnre. They suggested that a non-hydrostatic core or 
a subsurface ocean are the two possible interior models of Mimas that are 
consistent with their observations. It is evident that age can be calculated 
only if there exist measures of the number of observed craters. For example, 
in the case of Mimas, as there are no counts for craters with D < 3.5 km, 
the age-curve in Figs. [3] and |5] begins at this diameter. We calculate the 
age depending on D from Eq. |5l considering crater counts of the equatorial 
(E), south polar (SP) and intermediate latitudes regions (T) from Plescia 
and Boyce (1982), and also crater counts by KS09 from Cassini images of 
half of Mima’ surface. We plot the results in Fig. O Crater counts for the 
equatorial zone in Voyager images are given only for large diameters due to 
low resolution, and the age is nearly 4.5 Gyrs old. 

The south polar terrain of Voyager images has a similar age curve to 
Cassini images for D <12 km, but for larger craters, the south polar terrain 
seems to be younger, probably due to the lack of craters > 20 km found by 
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Plescia and Boyce (1982). Since on the one hand, Cassini connts inclnde 
almost half of Mimas’ snrface and have a better resolntion and on the other 
hand, SP area is small, it is possible that this drop of the SP age-cnrve is 
very local. For T-area and Cassini images, craters with D < 15 km are 
yonnger than 4.5 Gyrs old, and smaller craters might be even yonnger. For 
example, craters with D < 6 km are yonnger than ~ 3.5 Grys old. In Fig. 
m there are ~ 8000 craters per million km^ with D > 6 km that conld 
have been erased from the snrface in the age of the Solar System. Therefore, 
one possibility is that, effectively, small craters conld have been erased dne to 
particle depositions from the E-ring or even from a low-level geologic activity 
that shonld not be discarded on Mimas. However, since Mimas is heavily 
cratered, another possibility is that for D <20 km, crater density has reached 
a satnration eqnilibrinm. In fact, KS09 fonnd that the power-law index of the 
observed cnmnlative crater SFD is —1.548 for 4 km< H < 10 km and —2.12, 
for 10 km < H < 20 km, both indexes near to the satnration eqnilibrinm 
condition. The slope behavior of the age-cnrve is similar to Tethys and both 
cnrves are below the other two cnrves of the satellites plotted, which shows 
an erasnre process or satnration that acts on a wider range of crater sizes on 
those satellites. 

Becanse of the recently snggested potential for activity on Mimas, we 
cannot discard an erasnre process acconnting for the pancity of small craters 
snggested in onr model. However, the more conservative hypothesis would 
be saturation equilibrium for craters less than 20 km in diameter. 

3.3. Tethys 

Tethys is cold and airless, 1072 km in diameter, and very similar in size 
and aspect to Dione and Rhea. As it has a density of 0.9855'r/cm^ (Thomas, 
2010), very similar to liquid water, it is probably almost entirely composed 
of water ice and a small amount of rock. In fact, spectral observations of 
Tethys’ surface show no other component than H 2 O ice (Emery et ah, 2005). 
This satellite has two main features on its surface: one of them is a great 
impact crater known as Odysseus, whose diameter is ~ 450 km, and the other 
one is a great graben called Ithaca Chasma, whose length is ~ 2000 km. The 
interior of Odysseus seems to be younger than the rest of Tethys’ surface. In 
general, the surface of Tethys seems to be old and highly cratered, showing 
little geologic diversity. 

Royer and Hendrix (2014) present observations from the Ultraviolet Imag¬ 
ing Spectrograph subsystem of Cassini for Mimas, Tethys and Dione. They 
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found that Tethys and Dione have a leading hemisphere brighter than their 
trailing hemisphere at far-UV wavelengths, and Tethys shows a narrower op¬ 
position effect, reflecting a more porous surface. They attribute this structure 
to an intense bombardment by E-ring. 

There are crater connts on Tethys’ areas from Voyager images (Plescia 
and Boyce, 1982 and 1983) and from Cassini images (KS09). From those 
observed crater counts and our model, we calculate the age depending on D 
from Eq. |5] for those terrains. The two bands, along the west side of the 
terminator observed by Voyager and stndied by Plescia and Boyce (1982) 
have ages near 4.5 Gyrs, consistent with the high craterization observed. The 
other fonr zones analyzed by Plescia and Boyce (1983) have crater connts for 
D > 10 km. The heavily cratered terrain near 60° longitnde is 4.5 Gyrs old. 
The remaining three areas that include two plain nnits (GT and P) and the 
interior of Ithaca Ghasma graven (IG) have craters with D < 15 — 20 km that 
are yonnger than 4.5 Gyrs old, which is consistent with the age calculated 
by Plescia and Boyce (1985). Grater connts on Gassini images performed by 
KS09 were done for ~ 0.2 < H < 90 km and they have similar ages to plain 
nnits and IG for H > 10 km. Ages depending on D for GT, P, IG and KS09 
terrains are plotted in Fig. [6] and the heavily cratered terrains have been 
omitted for clarity. 

In Fig. [3l we plot the age-curve of Tethys only for KS09 crater connts, 
since it has a wider range of diameters and a better resolntion. It shows an 
age-curve similar to that of Mimas for the range of diameters plotted. Graters 
with H < 20 km could be young and, like the other satellites, smaller craters 
are younger than larger ones. KS09 and KSIO calculated the cumulative 
slopes of the SFD of craters and fonnd that for 0.2 km < H < 10 km the 
cumulative slope is —1.728, and for 10 km < H < 60 km is —2.2. Then, for 
small craters those distribntions could be close to the saturation eqnilibrinm, 
especially for 77 < 10 km. Graters greater than 77 ~ 20 km seem not to 
be affected by the erasing process or satnration, and thus the surface has 
preserved them, as can also be seen in Fig. 5 of DZ13. Therefore, those 
craters should be primordial. 

Like on Mimas, observations by Bnratti et al. (2011) do not show plume 
activity, and the npper limit of them are also one order of magnitnde less than 
the production on Enceladus. However, this does not rule out some possible 
activity. Since there appear to be small craters that would have been erased, 
there could, in fact, be some kind of geological activity (even plume activity 
at a low level, which has not been detected yet) that is acting on the satellite. 
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Figure 6: Age of the Ithaca Chasma (PB-IC) interior terrain, plain units (PB-CT) and 
(PB-P) from Voyager images and Cassini images (KS09) of Tethys with respect to D. 


and/or also bombardment of the E-ring flnx, as was snggested by Royer and 
Hendrix (2014). As can be seen in Fig. |2l this process should be able to 
erase, for example, ~ 14000 craters per 10® km? with H > 4 km or ~ 10® 
craters per 10® km^ with D >1 km in 4.5 Gyrs. 

However, as mentioned before, crater saturation in small crater sizes is 
the most conservative hypothesis for the paucity of observed counts. 

3 . 4 . Dione 

Dione has a diameter of 1122.8 km and a mean density of lA78gr/cm^. 
Its surface is uneven, as it presents heavily cratered terrains with craters 
greater than 100 km in diameter, as well as smooth plains (Smith et ah, 
1982; Plescia and Boyce, 1982). It shows a marked asymmetry between 
its leading and trailing hemispheres, with the leading side brighter than the 
trailing one. Royera and Hendrix (2014) analyzed the photometric properties 
of Dione revealing a more absorbent surface. They ascribe this to the fact 
that Dione is less intensely bombarded by E-ring grains and thus it should 
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have less fresh bright water-ice on its surface; and/or to that an exogenous 
darkening agent could be acting simultaneously on the trailing side. 

In addition, Buratti et ah (2011) analyzed spectral observations of Dione 
in search for plume activity and found that the upper limit of water vapor 
production is two orders of magnitude smaller than that of Enceladus, sug¬ 
gesting that this world is inert. However, plume activity cannot be absolutely 
ruled out because it can be under the detection limit. Further analysis by 
the Cassini magnetometer suggests that Dione has a tenuous atmosphere of 
transient origin (Simon et ah, 2011). Future observations by Cassini could 
possibly conhrm the existence of an atmosphere and/or some plume activity. 

Plescia and Boyce (1982) analyzed three terrain types observed by Voy¬ 
ager: a heavily cratered rough terrain (HC), an intermediate crater plain unit 
(IC) and a smooth plain unit (LC). They obtained cumulative size frequency 
distributions of craters for those terrains, which we used in our model to 
calculate surface ages. The HC area has crater counts for ~ 10 < D < 100 
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km and is 4.5 Gyrs old. Ages depending on D of IC and LC plain units are 
shown in Fig. [71 The IC unit has crater counts for 6 < D < 20 km and it 
seems that craters between 15 and 20 km are slightly younger than 4.5 Gyrs 
old. The LC terrain would be younger than 4.5 Gyrs old for all the diameters 
observed. This is consistent with the low crater density of LC and IC plains 
with respect to HC terrain found by Plescia and Boyce (1982) and also with 
the younger ages found by Plescia and Boyce (1985) for those IC and LC 
plains. 

On Dione there are two different surfaces recognized by KSIO from Cassini 
observations. The cratered plains (cp), which are heavily cratered, presenting 
records of primordial bombardment, and the smooth plains (sp) with less 
cratering, suggesting a younger surface. In Figs. [7] and El we plot the age 
with respect to D for both terrains on Dione. 

In these hgures we can see that the observed cp surface portion of Dione 
is generally old, but craters with D < 5 km could be young, and smaller 
craters would be younger than larger ones. This could be the result of a 
process in which craters have been gradually erased. This possible erasing 
process does not appear to affect craters greater than ~ 5 km. Although 
there is no data for D < 4 km for sp, the smooth terrain has a lower amount 
of craters than the cp surface, and the corresponding age curve is below the 
cp-age curve. For example, for craters with D > 4 km our model suggests 
that sp has an age of 3.8 Gyrs old while the cp is 4 Gyrs old. Indeed, since 
the age of sp for D > 8 km is near 4.5 Gyrs old (Fig. E]) and the number 
of calculated craters is greater than the observed one, this suggests a very 
slow erasing process for D > 5 km, which is not present in the case of cp. 
Potential erasing process/es on Dione would need to be able to erase in sp 
terrain, ~ 5500 craters per of D > 4, while in cp terrain there would 

be ~ 4500 craters per 10® km^ of D > 4 erased in 4.5 Gyrs. Comparing 
age estimations between Voyager and Cassini plains, it is evident that LC 
plains would be younger than smooth plains, at least in the range of diameter 
observed by both missions. It could be that LP plains were in fact young 
due to, for example, the internal activity which resurfaced portions of the 
surface of Dione, as was suggested by Plescia and Boyce (1985). However, 
Cassini images have a resolution of 1 to 2 orders of magnitude higher than 
Voyager images, so the comparison between age calculation with the data of 
both missions has to be taken with caution. 

On the other hand, like on the other satellites, KSIO calculated the cu¬ 
mulative crater SFD for cratered plains and obtained a cumulative slope of 
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b = —1.640, for 0.25km < D < 4 km, b = —1.166 for 4 km < H < 10 km, 
b = —2.31 for 10 km < H < 30 km, and —2.9 for 30km < D < 150 km. 
Therefore, the crater SFD of cp is similar to our model for craters greater 
than 30 km, but for smaller sizes the slope of the SFD is getting more flat¬ 
tened, similar to what occurs for lapetus, indicating saturation equilibrium 
for D < 10 — 30 km. For the smooth plains a similar behavior of the crater 
SFD has been found. Plescia and Boyce (1982) found a cumulative slope of 
b = —2 for the heavily cratered terrain, indicating saturation equilibrium, 
and a cumulative slope of b = —4 for the younger units, more similar to our 
model. 

Then, crater saturation could also be responsible for the paucity of small 
craters on Dione although there is morphological evidence that small craters 
in the sp have also been affected by endogenous resurfacing (Schenk and 
Moore, 2009). 

3.5. Rhea 

Rhea is the second largest moon of Saturn, with a diameter of 1529 km 
and a density of 1.237gr/crrR (Thomas 2010), a value somewhat greater 
than that of liquid water. Therefore, it is thought to be composed of a 
homogeneous mixture of ice and rock. Its surface of Rhea is densely cratered 
with abundant large heavily degraded craters. Rhea seems to be more heavily 
cratered than Dione and Tethys. An outstanding feature of this moon is a 
fresh ray crater of ~ 47.2 km of diameter, known as Inktomi. This crater is 
probably the youngest feature on Rhea’s surface, whose estimated age ranges 
from 8 to 280 Myrs old (Wagner et ah, 2008). Measurements by Cassini 
spacecraft have detected a tenuous atmosphere of oxygen and carbon dioxide, 
which appears to be sustained by chemical decomposition of the surface water 
ice under irradiation from Saturn’s magnetospheric plasma (Teolis et ah, 
2010). Stephan et al. (2012) analyzed the spectral characteristics of Rhea and 
concluded that the major process that affects the surface properties of Rhea 
is the interaction between the surface material and the space environment 
which includes the impacts of energetic particles from the magnetospheric 
plasma. This exogenous process could be responsible for the erosion of small 
craters on the surface of Rhea. 

Plescia and Boyce (1982) analyzed different areas of Rhea’s surface and 
found signihcantly different crater density. From our model, we found an age 
of 4.5 Gyrs old for the north polar region, this zone being densely cratered. 
Equatorial regions show subdued surface texture and absence of small craters. 
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Figure 8: Age of equatorial mantled (PB-M) and unmantled regions (PB-UM) from 
Voyager images and Cassini images (KS09) of Rhea with respect to D. 


suggesting that some areas have been mantled. Thereby, Plescia and Boyce 
(1982) separated them into mantled (M) and unmantled regions (UM). KS09 
obtained crater counts from Cassini images with a resolution of 0.01 — 1 
km/pixel for ~ 25% of Rhea’s surface, allowing them to reach very small 
diameters. We calculated ages for this Rhea’s surface and plotted them 
together with the equatorial regions analyzed by Plescia and Boyce (1982) 
in Fig. IHl We can see that the mantled regions are younger than the other 
terrains, that is, ~ 3.8 Gyrs old for H > 10 km, consistent with the age 
estimations by Plescia and Boyce (1985). However, although it is possible 
that this zone could be younger for smaller sizes, it cannot be guaranteed 
since there are no counts for craters smaller than 8 km in diameter for these 
zones. Then, in order to compare Rhea’s surface age with the other satellites, 
we plot our results for KS09 counts in Fig. [3l 

Rhea has the oldest surface with respect to the other satellites, with the 
exception of lapetus (Fig. [3]). However, craters with D < 4.5 km could be 
young, and the slope of the age curve is steep, which causes age to strongly 
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decrease in direct proportion to crater size. Thus, the possible process that 
might have erased small craters has a strong dependence on size. This agrees 
with the above mentioned process that gradually erodes the surface of Rhea, 
thus erasing small craters without burying the greatest ones. In Fig. [2l 
for example, there are ~ 40000 per 10® km^ craters with D > 2 that are 
erased from the surface of Rhea in the age of the Solar System. A peak near 
~ 0.5—km craters can be seen in the age curve. For D < 1 km, the number 
of observed craters on Rhea seems to have increased from the theoretical 
slope (see Fig. 7 in DZ13). Therefore, it is expected that our calculated age 
curve will change its behavior. As KS09 noted, for their observed impact 
crater distributions on the Saturnian satellites, D < 1—km craters are likely 
to be contaminated with secondary craters at some level, thus this “peak” in 
the age curve is probably due to this contamination on Rhea in the observed 
crater counts for small diameters. 

In relation to the observed cumulative crater SFD, KS09 and KSIO ob¬ 
tained cumulative slopes near —2 for all the diameter ranges. Then, although 
our model with S 2 = 3.5, that produces a crater SFD with a power-law index 
of ~ —3, hts better the number of observed craters, the power-law index 
of the distribution does not. Therefore, there exists the possibility that the 
surface of Rhea, is also saturated for the crater range observed. This would 
be consistent with the high density of craters observed and with the results 
of Squyres et ah (1997), who found that at least 25% of the craters on Rhea 
were destroyed by subsequent obliteration. 

3.6. Comparison with other age estimations 

As mentioned before, KS09 obtained the number of craters for the mid¬ 
sized icy satellites, and from previous cratering rate estimations by Zahnle et 
al. (2003), they calculated surface ages for each satellite and for some crater 
diameters too. Zahnle et al. (2003) analyzed two cases of a population of 
impactors. In case A, they inverted crater counts on the Galilean satellites to 
obtain the size distribution of impactors on Jupiter. In case B, they inverted 
crater counts on Triton from Voyager images to obtain the size distribution 
of impactors. In Table [21 we show both our age estimations for Mimas, 
Tethys, Dione and Rhea, for some diameters, and the estimations by KS09. 
The tendency of an increasing age with diameter that we have found is only 
shown by the results of case B (in general). Our results lie, in general, 
between case A and B, but for small craters they are more similar to case B. 
This could be due to the fact that the production of craters in the Neptune 
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Table 3: Ages obtained from our results and those calculated by KS09 all in Gyrs for 
various diameter ranges. 


Satellite 

D [km 

Age 

KS (A) 

KS (B) 

Mimas 

[7 

3.1 

4.4 

0.8 


10 4.2 

4.4 

1.3 


20 4.5 

3.3 

0.9 

Tethys 

1 

1.1 

4.6 

1.6 


r 

1.6 

4.6 

1.2 


[7 

3.5 

4.6 

1.7 


10 4.2 

4.5 

2.1 


20 4.4 

4.3 

2.7 

Dione-cp 

1 1.2 

4.6 

0.9 


2 

2.5 

4.6 

1.4 


[7 

4.4 

4.6 

2.6 


10 4.5 

4.6 

3.3 


20 4.5 

4.5 

3.1 

Dione-sp 


4.2 

4.6 

2.0 


10 4.5 

4.4 

2.0 


20 4.5 

3.3 

1.4 

Rhea 

1 1.6 

4.6 

1.7 


r 

3.4 

4.6 

2.2 


[7 

4.5 

4.6 

3.1 


10 4.5 

4.6 

3.7 


20 4.5 

4.5 

3.9 


System could be more directly related to SDOs just entering the Centaur 
zone, which is the case we are working on. However, the estimations by 
Zahnle et al. (2003) have an uncertainty factor of 4. 

According to our results, all surfaces analyzed up to now appear be old. 
However, we have noticed that for small sizes the observed number of craters 
is smaller than the calculated number from our model. The size limit of each 
satellite is presented in Table [H Although this size limit is different for each 
satellite, a common explanation for this paucity of small craters could be 
crater saturation. In particular, crater saturation starts in small sizes and 
would gradually incorporate larger crater sizes. 
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3.1. Enceladus 

Enceladus is a very interesting and beautiful moon of Saturn. It has a 
diameter of 504 km and a density of 1.609gr/cm^ (Thomas, 2010). Cassini- 
Huygens images reveal a very active surface with unique features, which 
invites one to investigate them. The discovery of the plumes of water vapor 
and small icy particles in the south polar region by Cassini Huygens in 2005 
(Porco et ah, 2006) showed us a very active world and at the same time, it 
was very important for the study of icy satellites. High resolution images of 
the south polar terrains from Cassini reveal tectonic features, a region almost 
entirely free of impact craters but with several ~ 130-km-long fractures called 
“tiger stripes”, probably warmed by internal tidally generated heat (Porco et 
ah, 2006). It was also found that this south polar region is 20 K hotter than 
expected from models, and a thermal emission of 3 to 7 GW was detected 
(Spencer et ah 2006). Gravity held studies by less et al. (2014) demonstrated 
that the structure of Enceladus is compatible with a differentiated body 
with, for example, a relatively low core density of ~ 2.4gr/cm^ and an H 2 O 
mantle of density of Igr/cm^ and thickness of 60 km. However, the local 
south polar heat huxes, gravity and topography, is more consistent with a 
regional sub sea of ~ 10 km thick in the south pole located beneath an ice 
crust ~ 30 — 40 km thick (less et ah, 2014). The surface of Enceladus is 
covered in almost pure water ice (Brown et ah, 2006). Newman et al. (2008) 
mapped the distribution of crystalline and amorphous ices on the surface 
of Enceladus from photometric and spectral analysis of data from Cassini 
Visual and Infrared Mapping Spectrometer. They found a mostly crystalline 
ice surface with a higher degree of crystallinity in the “tiger-stripe” cracks, 
and amorphous ice between these stripes. “Ice blocks” were also observed 
on Enceladus’ surface, like the regolites in asteroids and airless bodies of 
the Solar System (Porco et ah, 2006). However, their origin seems to be 
different. Martens et al. (2015) analyzed their distribution and origin in 
detail and found that, although they are found in other regions, they are more 
concentrated within the south polar terrain, where they must be produced 
mainly by tectonic deformation. 

All this suggests a very active and evolving surface. Moreover, Enceladus 
displays different types of terrain; in particular, as one moves away from the 
“tiger stripes”, more craters are found (KS09). This is due to the continuous 
resurfacing in the south polar zone. In fact, Porco et al. (2006) estimated 
ages within the south polar terrain possibly as young as 500000 years, or 
even younger. 



As in the previous section, we estimate the age of Enceladus’ surface 
comparing our calculations of the cratering production with the observed 
crater counts. There are crater counts on Enceladus from Cassini by KS09, 
but also previous counts from Voyager by Plescia and Boyce (1983) and 
Smith et ah (1982). Images from Voyager spacecrafts 1 and 2 showed a 
wide diversity of terrains on the surface of Enceladus. Different areas have 
been identihed, according to the populations of craters and other landforms 
(Fig. 21 of Smith et ah, 1982). Among these, the cratered terrains (ctl) and 
well preserved terrains (ct2) are distinguished. Another kind of terrain are 
the cratered plains (cp), whose principal landforms are bowl-shaped craters. 
In the equatorial region Voyager observed smooth plains spl and sp2. The 
grooves are the principal landforms in spl, which have few craters of 2 km 
to 5 km in diameter, while sp2 and ridged plains (rp) are almost free of 
craters. The number of observed craters in Voyager images were taken from 
Fig. 1 and Table 1 in Plescia and Boyce (1983) for ctl and ct2 (ECT) and 
cp (ECP), and from Fig. 23 in Smith et ah (1982) for spl, sp2 and rp. 

Kirchoff and Schenk (2009) analyzed Cassini images for crater counts. 
They analyzed different types of terrain separately according to variations in 
crater density and geological features (see their Fig. 4). Those regions are 
cratered plains (cp-eq and cp-mid) and the ridged regions rpl, rp2, rp3, rp4, 
rp5, and rp6. 

The regions analyzed by Voyager and Cassini seem to be different. How¬ 
ever, through an inspection of cartographic base maps of Enceladus in both 
papers (Smith et ah, 1982 and KS09), it seems that rp and sp2 and probably 
spl regions of Voyager are near the ridged plains (rpl-rp6) of KS09; ctl and 
ct2 regions of Voyager are near mid-latitude cratered plains of KS09, and cp 
regions of Voyager are near equatorial cratered plains of KS09. It is worth 
noting, however, that the Cassini data have a higher spatial resolution than 
that from Voyager. The Cassini mission has been able to show craters with 
D < 2 km. Also, KS09 has made a subdivision of ridged plains (rp), resulting 
in better details. 

From our model developed in DZ13, we have got the cumulative number 
of craters produced by current Centaurs. In order to visualize all the observed 
crater counts and to compare them with our model, those numbers are shown 

in Fig. [9l 

In general, the observed curves are below the calculated curve in our 
model. Therefore, the observed number of craters on the studied terrains 
would be lower than the current expected contribution of Centaurs with the 
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Figure 9: Cumulative number of craters calculated from model by DZ13 (s 2 = 3.5) and 
crater counts from Voyager (ECP, ECT, spl, sp2-rp) and from Cassini by KS09. 
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Figure 10: 


Age of the cp and rp6 surfaces of Enceladus with respect to D. 


exception of ECP and ECT, which present a greater number of observed 
craters than our model for D > 20 — 30 km respectively. 

We have calculated the age with respect to D from Eq. O and plotted 
the results in Figs. [10] and [11] for both observed crater counts (Voyager and 
Cassini) and regions. We have separated the results into two figures because 
of the wide age range. We can see that the observed areas of Enceladus 
are usually young, and that age increases as we move away from the “tiger 
stripe” cracks. Cratered plains in mid and high latitudes (cp-mid) are ~ 4.3 
Gyrs old, regions ctl and ct2 of Voyager are ~ 4.5 Gyrs old, but craters less 
than ~ 6 km could be younger. On the contrary, cratered plains in equatorial 
latitudes (cp-eq) are younger than 3 Gyrs old for the crater sizes counted by 
KS09. ECP from Voyager seems to be old, but craters less than ~ 20 km 
could be younger. 

All rp regions appear to be very young. The rpl region, which is very 
near the “tiger stripes”, a very active region, is the youngest region, and rp6 
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Figure 11: Age of the rp surfaces of Enceladus with respect to D. 
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is the oldest one. However, all sizes of craters in rp6 are young. In fact, 
there are craters with D <16 km which are younger than ~ 3.5 Gyrs old. In 
rp6, age increases with diameter like in the other mid-sized satellites. But in 
the rpl-rp5 zones, ages are similar for all the diameter ranges analyzed and 
they are all less than ~ 500 Myrs old. The very young region rpl is ~ 70 
Myrs old on average. Those ages of the ridged plains would be in agreement 
with estimations by Plescia and Boyce (1985) who estimated an age of no 
more than 100 Myrs. This behavior in the rpl-rp5 zones must be due to the 
plumes emanating from the south polar region and tectonic activity probably 
associated with “tiger stripes”, whose production mechanisms strongly erase 
those craters. 

In comparison with previous research, our ages are more similar to those 
calculated by KS09 for their case B (see Table 3). 

Moreover, Enceladus’ plumes are thought to be the source and mainte¬ 
nance of E-ring material, but a fraction of the plume particles returns to 
Enceladus hitting its surface (Kempf et ah, 2010). The particle deposition 
on Enceladus surface is more frequent in the “tiger stripes” than in other 
regions, as for example mid latitudes (Kempf et ah, 2010). Also, the largest 
water ice particles covering Enceladus’ surface (~ 0.2 mm) are concentrated 
in the “tiger stripes” zone and in general, the particle diameters are strongly 
correlated with the distance to this zone (Jaumann et ah, 2008). Therefore, 
the plume material returning to Enceladus could be an “exogenous” source 
causing the erasure of craters preferentially near the tiger stripes. 

On the other hand, the heat emanating from the south polar terrain, 
detected by Cassini’s Composite Infrared Spectrometer, suggests internal ge¬ 
ological activity (Spencer et ah, 2006). The heat source for the production of 
the plumes, resurfacing and heat flow in the south polar region of Enceladus, 
has not yet been sufficiently understood. The estimated radiogenic and tidal 
heating resulting from the orbital eccentricity of Enceladus are not sufficient 
to account for the energy emanating from the south polar zone (Porco et ah, 
2006; Meyer and Wisdom, 2007). Besides, the absence of internal activity 
on Mimas, which has an orbital eccentricity greater than Enceladus’, further 
constrain this as an unlikely source. Porco et ah (2006) analyzed other tidal 
heating mechanisms, such as mean motion resonances among the Saturnian 
moons, the 2:1 mean motion resonance with Dione being the most important, 
and secondary resonance heating, but those mechanisms are also inadequate. 
Squyres et ah (1983) suggested that the current rate of tidal heating might 
be sufficient to maintain a liquid sea under a crust of 10 km, but a previous 
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stronger heating would be required to initiate the melting. This previous 
process might have been that Enceladus had a larger eccentricity in the past 
(to obtain a greater tidal heating) or was captured into mean motion or 
secondary resonances with other satellites (Porco et ah, 2006; Meyer and 
Wisdom, 2008). Meyer and Wisdom (2008) demonstrated that Enceladus is 
probably near its equilibrium eccentricity and that the current equilibrium 
heating rate is not enough to trigger the melting. Therefore, any previous 
stronger heating requires some non-equilibrium behavior. 

4. Discussion and conclusions 

Based on both calculations of the cratering rate and the number of craters 
produced by Centaurs, and the comparison with observations, we have es¬ 
timated the surface age of each observed terrain on each mid-sized satellite 
of Saturn. The results are given for each satellite, but in general terms we 
can say that there are less observed small craters than calculated (except on 
lapetus), which result in young ages calculated from our model. This could 
be interpreted as either efficient endogenous or exogenous process(es) erasing 
small craters or crater saturation. 

We have also found that the smaller craters are being preferentially re¬ 
moved. This agrees with a gradual process of erosion, either from a geological 
or exogenous origin or by cratering itself. This erasing process does not affect 
large craters, and the size limit depends on the satellite. For craters greater 
than this size limit, our results imply that there are primordial craters, and/or 
that another main source of craters might exist. 

We have also found a correlation between the limit size from which small 
craters are erased (and then age is below 4.5 Gyrs old) and the distance to 
Saturn. That is, the smaller the distance from Saturn, the higher the size 
limit (see Tabled]). It is possible that there is not a general process to de¬ 
scribe this behavior since each satellite has its own peculiarities. However, a 
potential common process that seems to explain the paucity of small craters 
is crater saturation which could be connected to the correlation between dis¬ 
tances to Saturn and erasure size limit. Another possibility is that deposition 
of E-ring particles could be responsible for such a correlation. Since the E- 
ring has a density that generally decreases from Enceladus to Rhea (Horanyi 
et ah, 2008), the E-ring flux on the mid-sized satellites could be lower as we 
move away from Enceladus. 
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The estimated ages for the terrains analyzed on the mid-sized satellites, 
with the exception of Enceladus, are ~ 4.5 Gyrs old. Our ages are, within the 
ages calculated by KS09, but the fact that smaller craters are preferentially 
removed seems to be strong in our calculations. Almost all the Enceladus’ 
terrains could be young, with the exception of the mid and high latitude 
plains. The regions near the south polar terrain could be as young as 50 
Myrs old. On the contrary, all the surface of lapetus is old and it likely 
records a primordial source of craters, smaller craters being likely approach¬ 
ing saturation until D <2 km-craters, where saturation is complete. 

The ages calculated here are based on the current number of Centaur 
objects as contributors to craters. We have used the best known estimations 
from the literature for both this number and the size-frequency distribution. 
However, for small sizes, this SFD is uncertain and also unreachable with 
current direct observations. That is why the theoretical crater contribution, 
the comparison with observations and the determination of terrain ages in 
this paper should be carefully interpreted and not taken as absolute ages. 
They could serve, instead, to enhance both the understanding of the pro¬ 
cesses occurring on the satellites and the SFD for small sizes. Finally, the 
comparison of our age estimations with geological time scales on Saturnian 
satellites could help us determine what is still unknown today about those 
wonderful worlds. 
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